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Time-Frequency Resolved Coherent Anti-Stokes Raman Scattering (TFRCARS) 



5_j ■ was recently proposed as a means to implement quantum logic using the molecular 



ro-vibrational manifold as a quantum register [R. Zadoyan et al, Chem. Phys. 266, 
(2001) 323]. We give a concrete example of how this can be accomplished through 
an illustrative algorithm that solves the Deutsch-Jozsa problem. We use realistic 
molecular parameters to recognize that, as the problem size expands, shaped pulses 
must be tailored to maintain fidelity of the algorithm. 
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1 Introduction 



Logic operations for quantum computation (QC) between two or more quan- 
tum bits (qubits) have so far been experimentally demonstrated in liquid- 
state NMR [1], trapped ions [2], and cavity quantum electrodynamics [3]. 
Recently, Time-Frequency Resolved Coherent Anti-Stokes Raman Scattering 
(TFRCARS) in the molecular ro-vibronic Hilbert space, was proposed as an 
alternative approach [4] . Through measurements on room temperature iodine 
vapor, it was demonstrated that the elements sufficient for executing universal 
quantum logic, namely, the single qubit rotations and the two-qubit controlled- 
not gate, are naturally contained in TFRCARS. Here, we give an illustrative 
example of how such elements can be combined into a useful algorithm. We 
show how a TFRCARS interference experiment can use the molecular vibra- 
tional level structure to solve the Deutsch-Jozsa (DJ) problem [5]. The DJ 
problem has become a benchmark for demonstrations of algorithms on proto- 
types of quantum computers; having so far been used in NMR [6-11], linear 
optics [12], ro- vibrational molecular wavepackets in a pump-probe experiment 
[13], and excitons in semiconductor quantum dots [14]. 

TFRCARS is a four-wave mixing experiment in which a sequence of three 
non-collinear, short laser pulses (labeled P, S and P') are used to resonantly 
prepare and manipulate molecular ro-vibronic coherences. The P-pulse acts on 
a statistical state to prepare the first-order polarization P^^^ = c(^, 
where ^ represents the rotational, vibrational and electronic quantum indices 
(i) V, x) ■ The S and P' pulses then allow manipulation of the coherence through 
stimulated transitions to sequentially prepare second- and third-order coher- 
ences, consisting of complex superpositions determined by the radiation field. 
Transform-limited interrogation of the evolving third-order coherence can then 
be accomplished through the projective measurement of gated detection of 
the spectrally dispersed third-order polarization, P^^^ [4,15,16]. Energy and 
momentum conservation conditions, enforced through spectral and spatial fil- 
tering of P'^^^ , allow its detection without any interference from the incoherent 
background. The process can be mapped out as logic operations on the molec- 
ular register of eigenstates, by identifying the P pulse as initialize, S and 
P' pulses as process, and the gate pulse as readout. The approach has sev- 
eral attractive features with regard to QC. Since transitions involve electronic 
resonances, they can be carried out on fs time scales. A very large ratio of 
coherence to process time (> 10^) is afforded, even when we allow for de- 
coherence to be determined by the non-fundamental limit of inhomogeneous 
translational distribution (doppler broadening) at room temperature. Unlike 
NMR, which aims to manipulate effective pure states, in TFRCARS the sig- 
nal is only from the polarized sub-ensemble [4]. Very large superpositions, 
10^ — 10^ ro-vibrational eigenstates, can be manipulated with precision. Fi- 
nally, since the signal consists of a coherent radiation beam, information can 
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be transferred effectively and can be used to cascade operations. Small scale 
CARS-QC experiments, using several qubits, are currently being set up at 
Irvine. The limited number of qubits and process steps in the presently con- 
ceived algorithms, suggest applications in quantum communication [17], and 
simple quantum information processing tasks, such as quantum privacy ampli- 
fication and cryptography [18]. Even the simplest of propositions in quantum 
control of molecules [19] must contend with complications inherent in molec- 
ular state structures, such as vibrational level structures imposed by anhar- 
monicity of electronic potential energy surfaces. As such, principles of optimal 
control [20] will invariably be required in devising pulse shapes tailored for 
various implementations [21,22]. We highlight these considerations, by using 
realistic molecular parameters, intuitively obvious pulse shapes, and a clear 
computational task of solving the DJ problem. Finally, we note that in the 
short time after suggesting the mapping of 4- wave mixing onto QC [4], there 
have been demonstrations of the general approach by other groups as well 
[23,24]. 



2 Deutsch-Jozsa algorithm 



The DJ algorithm is one of the first examples of a quantum algorithm ex- 
hibiting a speedup compared to classical computers. The algorithm solves the 
following problem: "There is a finite set X — {x^}^^^ where N is even, and 
an unspecified function / : X — > {0, 1} that is promised to be either constant 
[f{xk) = or 1 V/c] or balanced [f{xk) = on exactly half the inputs]. De- 
cide with certainty whether / is constant or balanced in the least number of 
evaluations of /." The DJ algorithm has been extensively analyzed [25]. We 
now describe an alternative algorithm for solving the DJ problem, one that is 
optimized for an interference experiment. Let Ti be an A'^-dimensional Hilbert- 
space, where = \X\. Define a basis in Ti in one-to-one correspondence with 
the elements Xk^X: {\xi), \x2), \xn)}- 

(1) Prepare the equal coherent superposition = Y.k=i \^k)- 

(2) Define an operator Ff by its action on the basis elements as: Ff\xk) — 
(-l)/(^*)|a;i). Then: l^s) = F0i) = Ef=i (-l)^(^'=))lxi). 

(3) Make a projective measurement on \xi). 

The output (the signal) becomes: 

N 

5iv(/) = (Xl|*2) = E(-lP'=^ (1) 
k=l 
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Since f{xk) is either or 1, the terms in the sum are +1 or —1, respectively. It 
is clear that for / balanced, the sum contains an equal number of positive and 
negative terms. The signal in this case is zero. For constant /, all the terms 
have the same sign, so the absolute value of the signal becomes maximal. 
The signal therefore distinguishes the two types of functions and hence this 
modified algorithm performs the same task as the standard algorithm [25]. 
Furthermore, it too does so using just one /-evaluation.^ 



3 Scheme of proposed experiment 

We now propose a concrete experiment on iodine (I2) vapor to perform our 
alternative DJ algorithm. We illustrate the proposal through numerical simu- 
lations using accurately known molecular parameters. Vibrational levels rep- 
resent the basis states \xk) of Fig. 1 shows the iodine potentials for the 
ground (X) and excited (B) electronic states. We choose pulse widths and 
processing times < 1 ps, on which time scale the evolution of molecular ro- 
tations can be ignored, and therefore, we need only to be concerned with the 
vibrational levels on each potential. The CARS process is illustrated in Fig. 1 
through the relevant time-circuit diagram (Feynman diagram in inset). 

Initially the molecule is in its ground electronic state |X, 0)(0,X|, where the 
first (second) index denotes electronic (vibrational) quantum numbers. 

The P (pump) pulse induces the transition J2i\B,Wi){0,X\ ^ \X,0) {0, X\. 
A short enough P-pulse is chosen to prepare the desired number of ly-states 
with comparable amplitude. This step is the preparation of the equal coherent 
superposition state 

After a delay of an integer multiple of the vibrational period on the B state, the 
S (Stokes) pulse arrives to induce the transition J2j "^j) (0; ^1 ^ J2i\B, Wi) (0, X 
The delay is kept short (once or twice the vibrational period of the B state), 
so the S pulse acts on the same superposition as that prepared by the P-pulse. 
Otherwise, the S-pulse must be shaped to compensate for dispersion of the 
vibrational packet due to the anharmonicity of the B-state. The S pulse is 
designed to encode the function /. This is done by phase-shifting its spectral 
components (colors), corresponding to the {B,w) — > {X,4) transitions (or x^), 
using a phase- mask. We need only consider shifts and n shifts here, i.e., the 
mask element multiplies the spectral component Xk of the pulse by (— 1)-^*^^*=^ 
A mask with a specific sequence of phase factors ((— 1)-'^*^^^^ (— 1)./^(^jv)^ jg 
therefore in one-to-one correspondence with a function /. Note, the S pulse 

^ Our algorithm is similar to one described in [26], except that it does not use a 
qubit tensor-product structure. 
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provides a many-to-onc mapping of the B-state vibrations on the X-state vi- 
brations, and generates a broad coherent superposition \X, Vj){0, X\. How- 
ever, only the |X, 4)(0,X| coherence carries the processed phase information. 
This step corresponds to the apphcation of Ff, i.e., the computation of |^2)- 

The amphtude transferred to the |X, 4)(0,X| coherence is measured by ap- 
plying the P'(probe) pulse, and then collecting the time-integrated spectrally 
dispersed anti-Stokes radiation over a pre-selected transition. For the P' pulse 
which induces the J2k l-^j ^fc)(0, J2j 1-^) 'i^j)(0, transition, it is suffi- 
cient to select pulse characteristics that provide a one-to-one map of vibra- 
tional levels. Thus, for P' centered on the {X,4) — > {B,22) transition, as long 
as the spectral width of the pulse is narrower than the vibrational spacing near 
t; = 22 of 85 cm~^ (therefore a P' pulscwidth longer than 0.7 ps), we are guar- 
anteed the unique projection {X, 4) {B, 22). The desired signal is therefore 
obtained by dispersing the anti-Stokes radiation through a monochromator to 
detect the spectral amplitude on the {B, 22) — > (X, 0) transition. It is useful to 
note here that by limiting the detection window at the rotational band origin, 
where rotational recursions occur on a timescale t > 10 ps, it is possible to 
eliminate all potential interferences from rotational evolution. 

Let us summarize the scheme: 

1) Preparation: The P pulse prepares the equal coherent superposition. 

2) Computation: The phase-masked S pulse performs the Ff operation. 

3) Measurement: The |X, 4) {0,X \ coherence is measured as the time integrated 
spectrally-dispersed CARS signal on the (S, 22) — > {X, 0) transition, after ap- 
plication of the vibrationally selective probe pulse. This signal is proportional 
to the absolute value of the coherence |X, 4)(0,X|. The output characterizes 
/■ 



4 Pulse design 

We require an |X, 0) \B,w) \X,v) transition scheme with the following 
properties: a) For the general DJ algorithm w should assume values. Here 
we consider for definiteness the case of = 4, 8 (equivalent to two and three 
qubits). b) The Pranck-Condon (FC) overlap between X and B should be 
as large as possible and roughly constant for all {0,w) and {v,w) pairs, in 
order to maximize the efficiency of coherent transfer by the pulses, c) By 
choosing a large frequency shift between the P and S pulses it is guaranteed 
that the single time-circuit diagram of Fig. 1 describes the CARS process. 
The following values satisfy the above criteria: v — A; w — 20 — 23 ior N — A, 
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and w = 18 — 25 for = 8. The product of FC factors characterizes the 
overall strength of the 0) \B,w) 1-^)4) transitions. The product 
is largest and roughly constant around w = 22, which motivates our choice 
of the w values. In Fig. 2, we show the chosen pulse-shapes for P and S in 
the frequency domain for N — A. P is broad, covering to = 20 — 23. S is 
red-shifted, and phase-masked. Fig. 2 shows the effect of the (1,-1,1,-1) 
mask, corresponding to a balanced function /. Both the P pulse, and the S 
pulse before masking have durations of 50 fs. The P' pulse which must provide 
vibrational selectivity, is the same color as S but its duration is taken to be 
1 ps. 



5 Simulations 

We use Morse functions for the X and B electronic potentials with parameters 
=12550 cm.-\ =2.666 A, /3^=1.858A-i, L>f =4500 cm-\ rf =3.016 A, 
p^=1.850k-\ the energy difference between the minima of the potentials is 
15647 cm~^ (see Fig. 1). We calculate the vibrational eigenfunctions on both 
surfaces using the sinc-DVR method, and obtain the FC factors [27]. The 
time-evolution is explicitly integrated using the energy-representation. Time- 
ordering of the P and S pulses was not enforced: contributions from (P,S) and 
(S,P) sequences were added coherently. Only the Liouville pathway depicted 
in Fig. 1 was taken into account. No dccoherence mechanism was included 
(decoherence, which occurs with r > lO^^^s, can be neglected on the execution 
time-scale of 1-2 ps). After the application of the two pulses, the prepared 
\X, 4) {X, 0| vibrational coherence (which is an exphcit functional of the phase 
mask and the delay between the P and S pulses) preserves its magnitude 
A{f,T) = 4)(X, 0|| and oscillates with the {X,0) — {X,A) beat-frequency. 
By explicit calculation of the time-integrated CARS spectrum, using the per- 
turbative evaluation of P^^\t) as described in detail in Ref. [4], we verify that 
the experimentally observable {B, 22) — > {X, 0) spectral component is propor- 
tional to A{f, r). In what follows we simply consider A{f, r) as the signal. 



6 Results 

The signal ^(/, r) as a function of delay-time r between the P and S pulses 
was computed for all 2^ different phase-masks (i.e., for all possible different 

encoded / functions, including those that are neither constant nor balanced) . 
Fig. 3 shows the A^ = 4 signals as functions of the delay time between the P 
and S pulses for a constant and a balanced function. For the constant function 
(left) the signal is large at integer multiples of the (B) period, while the signal 
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is close to zero for the balanced function (right). This is therefore a clear 
experimental signature that distinguishes constant from balanced functions. 
Notice, however, that the signal from the balanced function is not exactly zero. 
In Fig. 4a, we show the signal A{f, 1) for all 16 different encoded / functions vs. 
S4{f) [Eq. (1)] (only 8 points are actually shown, since the signal is symmetric 
in ±S'4(/)). For balanced functions S4{f) = 0, so they should yield zero signal. 
For the two constant functions |S'4(/)| = 4, so they should give maximal 
signal. For the other functions |>S'4(/)| = 2, so they should give half the signal 
amplitude. In reality, as expected from Fig. 3, there is a spread along the signal 
axis, indicating that not all balanced functions yield exactly zero signal. Fig. 4b 
shows the results for = 8 case. Here the spread of signals corresponding to 
the same class of functions is significantly larger, thus reducing the fidelity 
of distinguishing constant from balanced functions. It is useful to introduce a 
formal fidelity measure in order to quantify the performance of the algorithm. 
Perhaps the simplest measure is the distinguishability D, which we define as 
the ratio of signal strengths 

^(/con) 



where the maximum is taken over all balanced signals. In the ideal case all 
^(/bai) = so D = 1. A second measure is the Pearson's correlation coeffi- 
cient r for a straight line interpolation between the signal A{f, r) and |S'iv(/)|, 
which also accounts for the spread in the other (neither constant nor balanced) 
functions. This measure is important in case one decides to use our algorithm 
to also distinguish these other functions. In Table 1 we present these fidelities 
for different values and for different delay times. The reasons leading to a 
degradation in fidelity are the following: I) Non-uniform FC factors, causing 
different \xk) to be transferred with slightly different amplitude by the P and 
S pulses. II) Spectral breadth of the laser-pulses, with the same consequence as 
I). These effects lead to non-uniformity in the prepared and projected coher- 
ent superpositions, thus misrepresenting and |^2) somewhat. Ill) Anhar- 
monicity of the {B) potential, leading to wavepacket dispersion on both the X 
and B potentials. In order to explore these effects, we performed simulations 
for the N — 8 case, where all FC factors were artificially set equal (a pulse 
achieving this effect can in principle be shaped by amplitude masking), and 
the P pulse was spectrally broadened, by reducing its duration to 10 fs. We 
refer to this as the "tailored" case. The results are presented in Fig. 4c, and 
the improvement in distinguishability is remarkable. This convincingly demon- 
strates the role of the non-uniform superposition in fidelity degradation. The 
solution, a spectrally broader P pulse together with tailoring the amplitudes 
of the P and S pulses to compensate for the non-uniformity of the FC factors, 
greatly improves the fidelity of the algorithm. To test the role of anharmonicity 
we also performed simulations at delays t — 0, 2tb- Obviously, at zero delay 
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anharmonicity cannot play a role. On the other hand anharmonicity should be 
more significant at r = 2tb than at r = r^. These expectations are confirmed 
in Table 1: a significant improvement in fidelity is seen when comparing the 
bare N — S results to the tailored results, but only for r — 0, Itb- In contrast, 
for T = 2tb the fidelity remains low. The degradation for large N values at 
large delay should be attributed to wavepacket dispersion due to anharmonic- 
ity. The simple solution is to apply the S pulse at r = or tb- The choice 
r = tb is actually preferable since interference from the time-reversed (S,P) 
sequence at r = 0, as well as higher order processes (6-wave mixing etc.) not 
accounted for in our simulations, lead to smaller signal strength. Optimization 
of the phases of the S and P pulses (chirping) is another option that can be 
used for compensating for anharmonicity effects [28]. 



7 Discussion and Conclusions 

The CARS-QC proposal [4] is a promising new implementation of quantum 
logic, with potential near-term implications for quantum communication. Our 
main purpose in this work was to test the feasibility of this proposal by study- 
ing in detail the implementation of a benchmark quantum algorithm. We chose 
the Deutsch-Jozsa (DJ) problem for its conceptual and technical simplicity, 
and devised a modified algorithm that only requires quantum interference, not 
entanglement. This allowed us to directly test aspects related to the prepara- 
tion of input vibrational superposition states and their dynamical evolution. 
To maintain fidelity, it is necessary to devise pulses that compensate for non- 
uniform Franck-Condon factors and counteract anharmonicity effects. These 
can be achieved through amplitude masking and frequency chirping using spa- 
tial light modulators [29], to optimally tailor the fields for a given choice of 
delay between pulses. Here, we have indicated that intuitively obvious pulses 
already allow a clear demonstration of how computational algorithms may be 
implemented in TFRCARS. The strategy in the laboratory demonstration of 
this, and similar algorithms involving rotation-vibration-electronic qubits, will 
be to rely on genetic algorithms to design pulses optimized for particular tasks 
[30-32]. 
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Fig. 1. X (ground) and B (excited) electronic potential energy surfaces of iodine, 
modelled as Morse potentials. Also shown: time-cicuit diagram of pump (P), Stokes 
(S), and probe (P') pulses, and relevant vibrational states. Inset: Corresponding 
Feynman diagram 
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Fig. 2. Spectrum of P and S pulses. S pulse for (1,-1,1-1) mask is shown. 
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Fig. 3. Signal A as function of delay time r for different masks. Time in units of 
(B) vibrational period, a) Constant function, b) Balanced function. 
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Fig. 4. Signal A as function of SnU)- a) N = 4,h) N = 8, c) "tailored" N 
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